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A method is developed for calculating effective sums of divergent series. This approach is a variant 
ON | of the self-similar approximation theory. The novelty here is in using an algebraic transformation 

with a power providing the maximal stability of the self-similar renormalization procedure. The 
latter is to be repeated as many times as it is necessary in order to convert into closed self-similar 
expressions all sums from the series considered. This multiple and complete renormalization is called 
self-similar bootstrap. The method is illustrated by several examples from statistical physics. 
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I. INTRODUCTION 



The most powerfull analytic tool for solving realistic, and therefore hard, problems in theoretical physics and applied 
mathematics is perturbation theory. However, perturbation series are notoriously badly behaved : in the majority of 
interesting cases they are divergent. A variety of mathematical techniques have been invented to assign a finite value 
to the sum of a divergent series. Such techniques are generally referred to as renormalization or summation methods. 
These methods are not only useful to theoretical and mathematical physicists but they are crucial because they 
provide a way to recover physical information from perturbative calculations. Probably the most common technique 
used to assign a meaningful value to a divergent series is the Pade' summation [1]. Using the latter one converts a 
formal power series to a continued fraction. Truncating this fraction at successive orders one obtains the rational 
functions called Pade' approximants. To reach a reasonable accuracy of such approximants, one usually needs to have 
i^i - tens of terms in a perturbative series. 

Recently, a method has been suggested [2-4] permitting to ascribe meaningful values to the limits of divergent 
, sequences by exploiting just a few terms of perturbative series. This approach, called self-similar approximation 
theory [2-4], is based on the following ideas. 

First, one has to incorporate into the considered sequence additional functions whose role is to renormalize the 
sequence making it convergent. These functions, because of their role, are named governing or control functions. The 
latter are to be defined from fixed-point conditions [2-4]. There are several ways of introducing such control functions. 
■ One and a natural way is to include them into initial approximation [5,6]. Fixed-point equations may be written in 
the form of the minimal-difference condition [5-7], or the minimal-sensitivity condition [8-14]. A condition close to the 
ON i latter type has also been used in the potential envelope method [15,16]. In two simple cases, of a zero-dimensional and 
one-dimensional anharmonic oscillator, the control functions were found analytically for an arbitrary perturbation 
order [17-20] by requiring the convergence of renormalized perturbation theory. Another way of introducing control 
functions is by using a scaled basis [21-24], in which the scaling parameters, playing the role of such functions, are given 
by analytical expressions with coefficients adjusted empirically from the convergence of a numerical iterative procedure. 
One can also define scaling parameters without using their analytical representation, directly from numerical search 
providing the convergence of an iterative method [25-29]. In the present paper we suggest a new way of introducing 
O ■ control functions, different from all variants mentioned. 

Another step in the self-similar approximation theory is to construct an approximation cascade whose trajectory 
is bijective to the approximation sequence considered [30-32]. In this paper we construct such a cascade not for the 
' sequence itself but for its transform. Of course, the idea of considering a transformed series instead of the initial one 
is not new. This is, e.g., the basis of the known Borel summation. Another example is the use of the Chebyshev 
transforms instead of straightforward power-series representation [33]. What is new in our approach is the use of a 
power-low algebraic transform with powers playing the role of control functions. 

Then the constructed approximation cascade is embedded into an approximation flow. Integrating the evolution 
equation of the flow, we obtain a self-similar approximation [2-4, 30-32]. 

In this paper we define the control functions from the principle of maximal stability of the approximation cascade 
trajectories. This is done by minimizing the absolute values of mapping multipliers, which is equivalent to a quasifixed 
point condition [30-32], since multipliers tend to zero when approaching a stable fixed point. 
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The plan of the paper is as follows. In Sec. 2 the algebraic transforms are introduced and the main steps of the 
self-similar approximation theory are sketched, not going into mathematical details which can be found in Refs.2-4, 
30-32, 34. In Sec. 3 the procedure of self-similar bootstrap is defined, consisting of multiple self-similar renormalization 
of all sums entering into given series. In Sec. 4 a particular case of the self-similar bootstrap is considered leading to 
a nice representation in the form of multiple exponentials. The following Sections illustrate the approach by various 
examples emphasizing the generality of the method that can be applied to problems of quite different nature. 



II. ALGEBRAIC TRANSFORMS 



Assume that we are interested in finding a function f(x) of a real variable x £ (— oo, +oo). Without loss of generality, 
the function f(x) may be considered to be real, since the case of a complex function can be always reduced to that 
of two real functions. Let perturbation theory give for the function f(x) approximations Pk(x) with k = 0, 1,2, ... 
enumerating the approximation number. The standard form of pu{x) is a series in powers, not necessarily integer, 
of x. The series can even include logarithms, since the latter, can always be presented, using the replica trick, as an 
expression containing a nonintegcr power. 

The algebraic transform is defined as 

P k (x, S )=x s p k (x), (1) 

with s real. This transform changes the powers of the series p k (x) changing by this the convergence properties of the 
latter. Effectively, the approximation order increases from k to k + s as a result of Eq.(l). The transform inverse to 
(1) is 

p k (x) = x- s P k (x,s). (2) 

To construct an approximation cascade, we proceed as follows. Define the expansion function x = x(f, s) by the 
equation 

Po(x,s) = f, (3) 

where Pq is the first available expression from Eq.(l). Substituting x(f,s) back into (1), we get 

y k (f,s) = P k (x(f,s),s). (4) 

The left -hand side of (4) represents a point of the approximation-cascade trajectory corresponding to approximation 
(1). The transformation inverse to (4) reads 

Pk(x,s) = y k (P (x,s),s). (5) 
The function (4) realizes the endomorphism 

y k (f,s) : R x R + -» R 

Consider the family {y k : k G Z + } as a dynamical system in discrete time. Since the trajectory of this dynamical 
system, according to (4) and (5), is bijective to the approximation sequence {P k }, this system was called [30-32] the 
approximation cascade. In order to simplify the consideration, let us pass from discrete time to continuous one. To 
this end, embed the approximation cascade into an approximation flow: 

{y k :keZ + }c{y(t,...):teR + }, 

which means that the trajectory {y(t, /, s)} of the flow has to pass through all points {y k (f, s)} of the cascade 
trajectory. 

The evolution equation 

d 

-g^y{t,f,s)=v(y(t,f,s),s) (6) 

for the approximation flow, where v{f, s) is the velocity field, can be integrated for an arbitrary time interval, say, 
from t = k — 1 to t = k* , which gives 
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Vk df 

-%— = k*-k + l; (7) 



here 

Vk=y(k,f,s), yl=y(k*,f,s). 

The upper limit in (7) corresponds, according to (5), to an approximation 

P:(x,s)=y(k*,P (x,s),s) . (8) 

The moment t — k* is chosen so that to reach the approximation (8) by the minimal number of steps. That is, we 
require that the right-hand side of (7) be minimal, 

t% = min(fc* - k + 1) . (9) 

Under condition (9), expression (8) is called the self-similar approximation. 

To find (8) explicitly, we need to concretize in (7) the velocity field v(f, s). This can be done by the Euler discretiza- 
tion of (6) yielding the finite difference 

Vk(f, s) = y k (f, s) - yk-i(f, s). (10) 
Thus, using (5), the evolution integral (7) can be written as 

P * df 

t% , (11) 



/ 



Vk(f,S 

where 

P k = P k (x,s), P* k =P* k {x,s) . 

When no additional restrictions are imposed, the minimal number of steps for reaching a quasifixed point is, 
evidently, one, 

abs min t* k = 1. (12) 

Additional restrictions can be of different types. For example, if the value of the sought function at some point xq is 
known, we may require that the found approximation would coincide at this point with the given exact value. Looking 
for an approximation in the class of functions with a prescribed symmetry is another way of imposing restrictions. In 
some cases the asymptotic behavior of the sought function at x — > and x — > oo may be available. Then requiring the 
correct asymptotic properties also plays the role of such additional constraints. In what follows, we shall concretize 
these variants by illustrating them with explicit examples. 

In this way, the sole quantity that is not yet defined is the parameter s of the transformation (1). Recall that 
our aim is to find an approximate fixed point of the cascade trajectory, a quasifixed point, which, by construction, 
represents the sought function. Therefore, the power s of the transform in (1) is to be chosen so that to force the 
trajectory of the approximation dynamical system to approach an attracting fixed point. Recall that s here is nothing 
but a kind of the control function, so, it is to be defined by a fixed-point condition. As is discussed in the Introduction, 
there are several forms of equations defining fixed points. Here we opt for a condition following from the analysis of 
fixed-point stability properties. 

Considering the mapping given by the approximation cascade, we may introduce the mapping multipliers 

d 

Vk(f,s) = -Tjj y k (f,s) . (13) 

This is related to the local Lyapunov exponent X k through the formula X k = £ In • Consequently, pL k ~ e Afcfe . If at 
increasing time, here at k — > oo, the trajectory approaches an attracting fixed point, then — » A < 0. This implies 
that the multiplier fi k — > 0, as k — > oo. Another quantity related to the multiplier (13) is the predictability time [35] 
which can be defined as r k w |A/s| _1 , or r k ~ |fc/ln \fj, k \ \ . This is the characteristic time during which the motion along 
the cascade trajectory effectively approaches a fixed point. When the latter is attractive, that is, when the limit of 
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the local Lyapunov exponent \ k , as k — > oo, tends to a negative value A < 0, then \i k tends to zero, and at the same 
time, larger absolute values |Afc| lead to smaller characteristic times r k . 

These properties show that closer we are to a fixed point, smaller is the absolute value of the multiplier (13). Hence, 
we may define the control function s as that providing the minimum of the multiplier. Instead of the multiplier (13), 
as a function of the variable /, it may be more convenient to pass to its image 

m k {x, s) = ^ k (F (x, s),s) (14) 

being a function of the variable x. Then the control function s = s k (x) is defined by the equation 

\m k (x,s k (x))\ = min\m k (x, s)\ . (15) 

s 

Because the minimization of the multiplier makes the trajectory more stable, we can call Eq.(15) the principle of 
maximal stability. And the so defined control function s k (x) can be termed the stabilizing control function, or, for 
brevity, the stabilizer. Notice, for comparison, that another definition of the fixed point would be to require the 
velocity (10) to be zero, which is exactly the minimal-difference condition [5,6]. After the stabilizer s k (x) is found 
from (15), we substitute it into (8) and, using the inverse transformation (2), we obtain the self-similar approximation 

f* k (x)=x- s ^P k *(x,s k (x)) (16) 

for the sought function. 

At the end of this section, let us note that the choice of control functions from fixed-point equations is rigorously 
justified when the fixed point is stable, that is if \/j, k \ < 1. When at some point k the trajectory becomes unstable, we 
have to stop the calculational procedure at the last stable point. Another possibility is to restructure the considered 
perturbation series. For instance, the instability of the procedure often happens when we are trying to construct 
self-similar approximation for a divergent function. Assume that we are dealing with such a function f(x) which 
diverges as x — > oo. Then, the self-similar renormalization procedure may become unstable at large x. To avoid the 
instability, we can either consider the function f~ 1 {x) or can rewrite the series in powers of l/x, treating the latter 
as a new expansion parameter. Usually, after this reexpansion the stability is restored. In the following, we shall 
illustrate this possibility by an example and will suggest a simple trick, giving the answer without the reexpansion, 
although being equivalent to the latter. 



III. SELF-SIMILAR BOOTSTRAP 



The procedure of calculating the self-similar approximations (16), starting from a perturbative series p k (x) is now 
completely defined. The renormalized quantity f k (x) must be, by construction involving the stability properties, a 
much better approximation to the sought function f(x) than the initial perturbative series p k {x). To improve the 
accuracy, we may repeat the self-similar renormalization applying it to other series that are left in (16). 

For illustrating this multiple self-similar renormalization, consider explicitly a perturbative series 

k 

Pk{x) = ^a n x n , a ^0, (17) 

n=0 

containing integer powers of x, although, as is mentioned above, the procedure works for arbitrary noninteger powers. 
Following Sec. 2, write the algebraic transform 

k 

P k (x,s) = Y / anX n+s (18) 

n=0 

of (17). As is seen, the transform (18) corresponds to an effectively higher perturbation order, k + s, as compared to 
the initial series (17), of order k. Eq.(3) for the expansion function x(f, s) now reads 

P {x,s)=a x s = f , (19) 

from where 

= (-) 1/s - (20) 
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The cascade-trajectory points in (4) become 



The velocity field (10) writes 



y k (f,s) = J2*n(^) n/S+1 - (21) 
n=0 a ° 



Vk (f,s) = a k (^) 1+k / s . (22) 

ClQ 



Calculating the evolution integral (11), with condition (12), we get the approximation (8) in the form 

P£(x,8)=P k - 1 (x,a)[l- ^j-A'-i (Mr /fe - (23) 



The stabilizer Sk(x) is to be found from the minimization of the multiplier 

k 



m k ( x ,s) = Y j ^{l + -)x n (24) 



from (14). Then we obtain the self similar approximation (16), 



r k (x)=p k . 1 (x)[l J^-^p^x)}-^ , (25) 



sa 

where s = s k (x). 

In this way, the self-similar renormalization led us from the initial perturbative series (17) to the renormalized 
approximation (25). The latter contains a perturbative series p k -i(x) of the lower order than the initial pk(x). This 
can be written as the relation 

f* k (x) = F k (x, Pk ^(x)) (26) 

showing that (25) depends on x and p k -i(x). We may repeat the procedure renormalizingpfc_i(a;) and getting f k _ 1 (x). 
With such a double renormalization, we come from (25) to 

fl*{x) = F k (x, fUix)) = F k (x,F k ^(x,p k ^(x))). (27) 

Here the doubly renormalized f k *(x) is expressed through Pk-2{x). Repeating the self-similar renormalization fc times, 
we obtain the fc-fold self-similar approximation 

fZ-*{x) = F k (x,F k - 1 (x,...{F 1 (x,a ))...)), (28) 

where we took into account that Po{x) = a . It may happen that (28) contains other power series. Then we may 
renormalize them as well. When all power series are renormalized, so that none of them is left unrenormalized, we 
have 

r k -*{x)^f k {x). (29) 

This complete procedure of the self-similar renormalization of all power series, resulting in an expression fk(x) 
containing none of them, will be called the self-similar bootstrap. 

IV. MULTIPLE EXPONENTIALS 

There are particular cases of the multiple self-similar renormalization yielding to a nice exponential representation 
of the resulting formulas. Here we consider one such a sufficient condition when all coefficients at perturbative powers 
are positive, a n > 0. Then the minimum of the multiplier (24) realizes at s — > oo. Taking this limit in (25) gives 

r k (x)= Pk ^(x)eM-x k )- (30) 
ao 
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Repeating the renormalization in line with (26) and (27), we get 

fk*{x) = Pk- 2 (x) cxp{ — (a k ^ 1 x k - 1 + a k x k )}. (31) 

Continuing so on, we obtain the fc-fold approximation (28) in the form 

fk"*( x ) = a o CX P{ — {aix + a 2 x 2 + ... + a k x k )}. (32) 

As we see, the fc-times repeated self-similar renormalization does not deliver us from the power series. Really, the 
fc-star approximation (32) is expressed through a part of the initial perturbation series (17), namely, through 

k 

Pk(x) -a = ^2a n x n . 

n=l 

With the notation 

k 

Pk(x) = ^2a' n x n , (33) 

n=0 

in which a n = a n+1 , n = 0, 1; 2, ...k, we may rewrite (32) as 

X i 

ft"*{x) = a exp{— p fc _i(x)}. (34) 
The power series p fc _ 1 (x) can be renormalized in our standard way giving the corresponding self-similar approximation 

f'k-i(x) = aoexp{^-p^'_ 2 (x)}, (35) 

in which 

k 

p'k( x ) = a « a; ™' a « = a «+ 2 • ( 36 ) 

n=0 

With this renormalization in mind, we transform (34) into 

ff*(x)=a cM-.f'k-i(x)}- (37) 
a 



Combining (35) and (37), we have 



fk*-*{x) = a exp{— aiexp{— p'^ 2 (x)}}. (38) 
a ai 

Converting k times all power series in the exponentials, with the use of the notation 

b = a , b k = — , k = 1,2,..., (39) 
dk-i 

we obtain the bootstrap self-similar approximation 

fk(x) = b eyLp(bixexp(b 2 xexp(...bk-ixe?ip(bkx)))...), (40) 



as is discussed in (29). 

In the case of small x — > the expression (40) yields 



fk(x) ~ c + c\x + c 2 x 2 + c 3 x 3 (41) 



G 



with the coefficients 



1 11 

c = b , c 1 =b b 1 , c 2 = b bi(b 2 + -h), c 3 = b bi(b 2 b 3 + -b\ + b x b 2 + -b\) 

i AO 



Substituting here (39), we have 

a\ a 2 ciia 2 a? 
c = a , ci=ai, c 2 = a 2 + - — , c 3 = a 3 + 1 h — ^ . 

This shows that the asymptotic, as x — > 0, behavior of (41) and (17) coincides up to the linear terms. For the 
higher-order terms, a n ^ b n for n > 2. Such a renormalization of the higher order expansion coefficients is typical 
of the self-similar approximation theory [2-4]. This renormalization allows to extend the region of applicability of 
self-similar approximations, with respect to a variable x, as compared to the initial perturbative series. 

It is worth mentioning that the multiple, or continued, exponentials of type (40) have been studied in mathematical 
literature for already more than two centuries, since Euler; and a number of references can be found in Refs. [36,37]. 
We have derived the form (40) following the multiple self-similar renormalization for the power series (17) with positive 
coefficients. Of course, for other sets of coefficients in (17) the final bootstrap approximation will not necessarily take 
the form of a multiple exponential, as in (40), but will be a mixture of exponentials and fractions, each expression 
being conditioned by the principle of maximal stability (15). 

Another way of getting a multiple exponential could be as follows. Consider a sequence {ipk(x)} of the terms 

ipi(x) = 6 cxp(6ix), 
(f 2 (x) = &o cxp(bix cxp(b 2 x)) , 

and so on, where the coefficients bk are given by (39). At the fc-step of the sequence {ipk(x)} we come to (40). However, 
this way is justified if the approximation cascade corresponding to the sequence {ipk(x)} possesses a stable trajectory 
leading to a quasifixed point representing (40). This approximation cascade can be constructed in the standard way 
[32,34]. From the equation <pi(x) — f we find x = x(f) which is 

s(/) = ^ln(-f). 

Then we define the approximation cascade as is described in Sec. 2. The corresponding trajectory {zk(f)}, consisting 
of the terms 

z k (f) = <pk(x(f)) 

is bijective to the approximation sequence {ipk(x)}. 

The stability of the trajectory {zk(f)} is checked by calculating the multipliers 



^iM) = [^^] f ^,un 142) 



and analyzing them with respect to the stability condition |Mfc(x)| < 1. 

At the end of this section let us remark that we have used here the word "bootstrap" in the generally accepted 
meaning, as a kind of a completely self-consistent procedure permitting to construct an explicit solution to a compli- 
cated problem. This term in the close meaning was used, for example, in constructing a self-consistent distribution 
over particle masses in high-energy physics [38] , and also in constructing the ^-matrices for two-dimensional conformal 
field theories [39]. 

V. ZERO-DIMENSIONAL ANALOGS OF FIELD THEORIES 

A. Non-Degenerate Vacuum 

Consider the thermodynamic potential of a zero-dimensional anharmonic model (see e.g [32]) represented by the 
integral 
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1 f°° 

J(g) = -= \ exp(-.T 2 - gx 4 )dx, (43) 

with the integrand possessing a single "vacuum" state, located at the point x = . The expansion of this integral in 
powers of the coupling parameter g, around the vacuum state leads to divergent series, 

J(g) ~a + bg + eg 2 + dg 3 + hg A + (44) 

3 105 , 3465 , 675675 

a — 1, b = , c — , a = , h = . 

4' 32 ' 128 ' 2048 

We apply here the self-similar bootstrap renormalization, guided by desire to perform as many renormalization 
steps, as possible. Write down the following set of approximations to the quantity J(g), analogous to the general form 
(17): 

Jo(g) = a, 

Ji(g) = a + h a, 

J 2 (g)=a + bg + cg 2 (45) 

M9) =a + bg + cg 2 +dg 3 , 

J 4 {g) = a + bg + eg 2 + dg 3 + hg A , 
together with the following local multipliers, that can be found from the general representation (24): 

( \ 1 , bl+s 
mi(g,s) = 1 H g , 

a s 

i \ / \ c 2 + s 2 

m 2 (g,s) = mi g,s H g , 

a s 

m 3 (g 1 s)=m 2 (g 1 s) + ^'^-^g 3 , (46) 

m^g, s) = m 3 (g, s) H g . 

a s 

The analysis of (46) shows, that in the last three cases the most stable trajectories are realized at s — > oo, and that 
in the first case, s — > oo also corresponds to a stable trajectory. Therefore the starting four steps of the self-similar 
bootstrap renormalization can be safely performed in the exponential form, leading to the intermediate renormalized 
quantity 

J 4 ****(.9) = aexp{^(bg + eg 2 + dg 3 + hg 4 )}. (47) 

Write down a set of approximations to the quantity J'(g) — bg + eg 2 + dg 3 + /ig 4 , appearing in the exponential of this 
expression: 

j'i(g) = bg, 
j' 2 (g) = bg + eg 2 , 
J'M = bg + eg 2 + dg 3 , 
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j'M = bg + cg 2 + dg 3 + hg\ 



with the following local multipliers: 

// \ c2 + s 

m 2 {g,s) = 1 + iy+Ii 9, 

m' 3 {g,s) = m' 2 {g,s) + ^^-g 2 , (48) 
o 1 + s 

ii \ ii \ hA + s , 
m 4 {g,s) = m 3 (g,s) + t— 
o 1 + s 

The analysis of (48) leads us to the conclusion that the exponential renormalization is optimal at every step and, 
following to the standard prescriptions of Sec. 4, we transform (47) into 

Jt***(g) = acxp{-g exp[-J-(c ff + dg 2 + hg 3 )}} . (49) 
a b 

Our routine procedure requires to renormalize now the quantity J (g) = eg + dg 2 + hg 3 , using the following approxi- 
mations: 

■h (flO = c 9, 

J 2 (a) = c 9 + d g 2 , 

J 3 (g) = cg + dg 2 + hg 3 , 



and analyzing the following multipliers: 



hi v d2 + s 

m 2 {g,s) = 1 + 



m 3 {g,s) = m 2 (g,s) + . 

We conclude that the most stable trajectory corresponds, in both cases, to the exponential summation, leading to the 
intermediate formula 

Jr*{g) = aexp{-gexp[^ exp[-(dg + hg 2 )}}}. (50) 
a b c 

The last step of the procedure, applied to the quantity j "(g) = dg + hg 2 , with the approximations set: 

J\ (g) = dg, 

J 2 (9) = dg + hg 2 , 

and with the multiplier 

mi \ h2 + s 

1712 (ff ' s) = 1 + dIT^< 

again should be performed with s — > oo, and the bootstrap program is completed: 

J A ( g ) = aexp {^ ff exp[^exp[^ 5 [exp(^ 3 )]]]}. (51) 
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Similar expressions follow when less terms from the initial expansion are taken into account: 

Mg) = acxp{^exp[^exp(^.g)]}, 

— be 
Mg) = aexp{-gexp[-g]}, 

Ji(g) = aexp{^y}. 

We had pointed out already, that the last expression corresponds to the stable, but not to an optimal, i.e. the most 
stable, trajectory. Analyzing m\ (g, s), we obtain the optimally rcnormalized expression: 

M(g) = a[i V^r s(9) , S (g) = -^-. (52) 

a x s(gj a + bg 

At the point g = 1, the following numbers are generated by the sequence Jj, i = 2, 3...: 

Ml) = 0.472 (JTo = 0.512), T 2 (l) = 0.991, = 0.473, = 0.991. 

We observe the two subsequences, with odd and even numbers, with values practically unchanged within each sub- 
sequence, probably embracing the correct result from below and from above, respectively. We can suspect that the 
corresponding sequence Ji possesses the two competing unstable quasiffxed points, i.e behaves chaotically and, in 
such situation, it is appropriate to use the self-similar approximation smoothed by the Cesaro averaging procedure [32], 
i.e. in our case, simply to take the average over the two neighboring members of each subsequence. 

This conjecture, is supported by the analysis of the corresponding sequence of the multipliers for the sequence Ji , 
as discussed in Sec. 4. From the initial approximation 

Mg) = f, 

one can find the expansion function 

b y/ira 

and, after the routine transformations, the following expression for the multipliers (42) can be obtained: 

Mi(g) = 1, 

M 2 (g) = $ 2 (g)^2{g) 

Mais) = *s(ff)*3($), 

M 4 (g) = $ 4 (<?)*4(<7), 

where 



$ 2(5 ) = ^Mexp(^), 
a b 



= -^exp(-<7exp(- ff )), 
a o c 



* / x Mg) , c ,d ,h w\ 
®i(9) = — — exp(-5exp(- 5 exp(- 5 ))), 
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and 



*2(0) = b + cg, 

*3(ff) = & + cgexp(^) + dg 2 cxp(^g), 

* 4 (fl) = b + c 5 exp(|gexp(|flr)) + d ff 2 exp(§ 5 ) exp(| 5 exp(§. 9 )) + 
+h 3 3 exp(f 5 )exp(^exp(^)). 

The following values are obtained at g = 1: 

Mi = 1, M 2 = -0.089, M 3 = 1.008, M 4 = -0.089, 

supporting our initial guess, that the approximation cascade behaves chaotically. After the Cesaro averaging, the 

sought value at g = 1 equals to, say, • /3 ( 1 )+ j4 ( 1 ) = 0.731, deviating from the exact value 0.772 with the percentage 
error —5.228%, an acceptable accuracy if to remember that the initial expansion (44) gives the percentage error 
~10 4 %. With the optimized J\ Q , we obtain even better estimate, 0.752, for the sought value, with the percentage 
error equal to —2.668%. 

B. Double-Degenerate Vacuum 

Consider the integral 

/oo 
exp(x 2 — gx A )dx, (53) 
-oo 

representing zero-dimensional field-theory with the integrand possessing the two maxima, located at the points 



w 



where g plays the role of a coupling constant. We intend to estimate this integral in the region of intermediate 
couplings g ~ 1. It was pointed out in [40], that any conventional expansion, in powers of g, or g^ 1 , is not sufficient, 
since it does not take into account the existence of those degenerate maxima, corresponding to the double-degenerate 
"vacuum". 

Within the framework of D— dimensional field theories, the existence of degenerate vacuum is taken into account, 
e.g., by means of the zero-energy instanton-anti-instanton solutions, and all further corrections to observables come 
from the excitations above the instanton-anti-instanton background and from interaction of all those quasiparticles. 
In our case we take into account the double-degenerate vacuum by means of the shift: 

x = x — A, 

then expand the integral around the two saddle points and apply the self-similar renormalization to the resulting 
asymptotic expansion in powers of a small parameter g 1 / 2 , continuing it to the region of g <~ 1. Following these 
prescriptions, represent the integrand in the vicinity of one of the saddle points in the form: 

exp(x 2 - gx A ) w cxp(i-) exp(-2A 2 ) CX p(A(g)X 3 ) « cxp(^) exp(-2A 2 )(l + A(g)X 3 + ...), 

and perform the integration, so that 

1(g) ^2eM^)(a + bg 1/2 + ■■■), a = ^2^l\ b = 2^' 2 . (54) 



Applying the self-similar renormalization, we readily obtain 



1 bq 1 / 2 

r( 5 ) = 2aexp(-)exp(^— ). (55) 
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Despite the absence of dynamics in the zero-dimensional case, I*(g) consists of two factors of different nature; one 
of them is non-analytic in the coupling constant, resembling the well-known "instanton" term within the framework 
of nontrivial D— dimensional field theories, the second is analytic in g 1 ' 2 and resembles a contribution from the 
excitations above the instanton-anti-instanton background. 

The percentage error for the renormalized calculated with respect to the exact 1(1) = 2.762, is 2.462%, and 

considerable improvement is achieved compared to the percentage error of the perturbative expansion (54), equal 
to-8.834%. 



VI. STRONG COUPLING REGIME 
A. Zero-Dimensional Case. 



Let us apply to (43) the so-called "strong-coupling" expansion, in powers of 1/g, with a quartic term of the integrand 
taken as a zero approximation, representing the integrand of (43) as follows 

x 4 

exp(-x 2 - gx 4 ) w exp(-.ga; 4 )(l - x 2 + — + ...). 

After integration we obtain the expansion in inverse powers of the coupling constant 

7/ n -i/4 i _3/4 -5/4 1-813 , -0.612 0.227 tro . 
J(g) w ag 1 ' 4 + bg 3 ' 4 + eg 5/4 + a = b = — , c = (56) 

y/TT V7T V 71 " 

We write down the following consecutive approximations to the quantity J(z), where z = g^ 1 ' 4 , 

J\(z) = az, 

J 2 (^) = az + bz 3 , 

J 3 (z) = az + bz 3 + cz 5 ; 

and the multiplier m2(z 1 s) = 1 + f f^pfz 2 reaches its minimal zero value at s = 0.019, being much smaller than the 
minimal value of the corresponding multiplier 7713(2:, s) = 7712(2:, s) + f f^pf^ 4 - Therefore, the renormalized quantity 
J|(z), will correspond to the more stable trajectory, than J|(z). Following the standard prescriptions of Sec. 3 we 
obtain 

r 2 {x) = ax{l —^ x 2 } -l^) + W, s(x) = -i±^ . 

a[l + s(x)\ a + bx^ 

The percentage error for the renormalized quantity J| (1), calculated with respect to the exact J(l) = 0.772, is 2.266%, 
and a considerable improvement is reached compared to the percentage error of the perturbative expansion (53) with 
only starting two terms taken into account, equal to— 12.208%. 
Represent (43) in a little bit different form: 

J(g) = g-^(a + bg- 1 ' 2 + eg' 1 + ...) = g-^a + 1(g)], (57) 

and write down the following set of approximations to </(<?) , using the variable y = g^ 1 ' 2 : 

Jiiy) = fry, 
My) = b y + c y 2 - 

The multiplier 7772(7/, s) = 1 + 1 1±| y, is minimal at s = 0, and | 7772(2/, 0) l< 1- The renormalized quantity J2* (y) can 
be readily written down: 



Recalculating now J* (g), we obtain J*(l) = 0.771, with the percentage error —0.13%, much better than the percent- 
age error 4.386%, given by the initial expansion (57). Even at small g — 0.21, the percentage error given by the 
renormalized expression, remains less than 1%; at the same time, the percentage error given by the initial expression 
reaches 43.538%. 
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B. One-Dimensional Case. 



Consider the dimensionless ground state energy e(g) of the celebrated quantum one-dimensional quartic anharmonic 
oscillator, closely connected to the so-called ip A model in the quantum field theory (see e.g. [3]). Here g stands for 
the dimensionless coupling constant, expressed through the parameters entering the Hamiltonian of the system by 
the known relation (see e.g [3,41]). The asymptotic expansion for e{g) in the strong coupling limit, corresponding to 
g — > oo, is known (see e.g. [41]) in the following form: 

e{g) ^ ag 1/3 + bg~ 1/3 + cg~\ a = 0.667986, b = 0.14367, c = -0.0088. (58) 

Let us, using the experience gained while considering the strong coupling limit of the zero-dimensional field theory, 
renormalize the last two terms of the expansion (58). Using the notation y = g^ 1 ^ 3 , we write down the following set 
of approximations for the quantity e — e — ag 1 / 3 : 

ei(y) = by, 



e2 (y) =by + cy 3 , 

with the multiplier m,2(y, s) = 1 + ff^fy 2 , possessing minimal value at s = 0, when g y 0.1. The renormalized 
expression can be obtained following the standard prescriptions of Sec. III. Returning to the initial variable we obtain 

h 3 / 2 

e " (s)=as "' + 7m^Tc (59) 

An accuracy given by e*(g), can be elucidated by comparison with the "exact" numerical results (see e.g. [41]). At 
g = 0.3, the percentage error, given by (59), is equal to —0.099%, at g = 1 it is —0.022%, and at g — > 200 the 
percentage error tends to zero. To our knowledge, these results are better than those obtained by other analytical 
methods. On the other hand, at small g = 0.001, an accuracy of the formula (59) is by far inferior, compared to many 
other analytical methods. The reason can be understood if notice, that e*(0) = 0.41048, strongly deviating from the 
known value 1/2, but being much better than the infinite value predicted by the initial expansion (58). We conclude, 
remarking, that using the effective time t* as an optimization parameter, determined from the condition e*(0) = 1/2, 
one can achieve better agreement of the renormalized formulae with the exact results in the region of small coupling 
constants. For the goal being pursued in the present paper, it is fairly enough to limit the discussion by formulae 
(59), designed for the intermediate and strong-coupling regimes. 



VII. EQUATION OF STATE 

A. System of Hard Spheres. 

We demonstrate below, how the self-similar bootstrap can be applied in the theory of equations of state for simple 
liquids. Consider the model system of hard spheres with the diameter d [42,43], where the empirical equation of state, 
connecting pressure p, temperature T, the number density n and the reduced density p = irnd 2 /6 7 is known: 



nkT (1 - p) 3 

The equation of state (60) agrees very well with the molecular dynamics results [42] . On the other hand, the theoretical 
virial formula according to Percus-Yevick [41,42], is given as follows: 

p 1 + p + p 2 -3p 3 

(61) 



nkT (1 - p) 3 

These two expressions almost coincide at low densities, e.g at p = 0.1, the percentage error of Eq.(61), as compared 
with (60), equals —0.18% ; while for the intermediate and high densities the agreement becomes very poor, e.g. at 
p = 0.5, the percentage error is —15.385%, and at p = 0.8 it equals —53.112% . 
Consider the regular part of (61), defined as r : 
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r = 1 + p + p 2 - 3p 3 , 



(62) 



as an asymptotic, low-density expansion for the "true" regular part r(p), and try to continue the expression (62) 
from the region of p <C 1, to the region of p < 1. It seems reasonable to use for renormalization only the last three 
terms from (62), since the constant term describes the ideal gas behavior and we are interested in the region of high 
densities. Let us write down the following consecutive approximations to the quantity r = r — 1 : 



n 



ri = p + p 2 , 



n = p + p 2 -zp A 

The multipliers are 

2 + s 



m 2 (p,s) = 1 + p 



1 + 



3 + s 

m 3 (p, s) = TO 2 (p, s) - 3p 2 — — , 

i- \ s 

m 'i(p, s ) = 1 - ^fr - ^- 

It is admissible to apply here the self-similar bootstrap renormalization in the form of the continued exponentials, 
since at every step of the procedure the exponential summation is performed along the stable trajectory. Following 
the standard prescriptions of Sec. 4, we obtain 

r{p) = pcxp(pexp(-3p)). (63) 

The multiplier M(p), corresponding to (63), is given by the expression 

M(p) = exp[pexp(-3p)] exp(-4p)(l - 3p), 

and is very small at p > 1/3, e.g. M(0.8) = —0.061, signaling the robust stability of the sequence of the continued 
exponentials (63). Recalculating 



_p_ r(p) + 1 

nkT (1 - pf ' 1 ' 

and comparing it to the empirical formula (60), we obtain that at p = 0.1, the percentage error equals —0.118% ; at 
p = 0.5, the percentage error is —4.061%, and at p = 0.8 it equals —3.516%. 

We see that the equation of state (64), obtained from the bootstrap self-similar renormalization, much better and 
more uniformly agrees with the computer experiment, than the initial virial expansion (61), over the entire range of 
densities. The agreement drastically, by 17 times, improves in the region of high densities. 



B. System of Hard Hexagons 

The model of "hard hexagons" represents a simple two-dimensional model of impenetrable molecules on the tri- 
angular lattice. The model allows an exact solution [44], and the phase transition from the liquid phase existing 
above the critical value of the so-called activity z c — 11.0917..., to the solid phase, existing below z c , is well studied. 
The equation of state, describing the dependence of the order parameter R on the activity-related parameter, can be 
written down in quite a complicated and not very convenient form [44] . On the other hand, the critical value of density 
p c at the point of phase transition is known too, and equals 0.27639. ..[44]. Independently, the high-density expansions 
of the order parameter in powers of the inverse activity z' = 1/z, or in powers of the high-density variable p' — 1 — 3p, 
were obtained [45]. Their quality is considered as very high, since the critical parameters could be determined from 
them with extremely high accuracy, using the Pade approximants in conjunction with some extrapolation methods 
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[45]. We present below simple expressions for the equation of state, obtained as a continuation of the high-density 
expansions of the order parameter up to the point of phase transition. 

The expansion of the order parameter in powers of p' is given as follows [45] : 

R=l- 3(p') 2 - Hp') 3 - 36(p') 4 - 159(p') 5 - (65) 

Let us apply to (65) the bootstrap self-similar renormalization based on the exponential summation at every step, 
and leading to the equation of state for the system of hard hexagons in the form of the continued exponentials: 

~ 53 

R(p') = exp(-p'cxp(3p'cxp(3 /9 'exp(4 )0 'exp(— p' ))))). (66) 

The function R(p') approaches zero at p'* — 0.170005(±1), corresponding to the p* = 0.276665, deviating from the 
exact value by 0.1%. Thus, the renormalized equation of state (66), agrees with the initial expansion in the region of 
p' <C 1 by design, and predicts the point of phase transition with very high accuracy. On the other hand, the form of 
continued exponential, can be justified aposteriori, analyzing the multipliers (42), where it can be shown, after some 
lengthy, but routine calculations, that M 5 (p') — ► 0, in the region of p' w p'*, i.e. the stability condition is satisfied 
along the trajectory, described by the sequence of approximations corresponding to (66), in the vicinity of the critical 
point. Similarly, using the known expansion of R up to the fifth order terms in z', the corresponding equation of state 
can be obtained. In this case we found the critical z'* ~ 12.1803(±1), deviating from the exact value by 9.8%. 



C. Polymer Coil 

The expansion factor a of the polymer chain, within the framework of a standard " beads-on-string" model, is 
conveniently represented as a function a 2 =a 2 (z) of the parameter z = 2(^) 3 / 2 N 1 / 2 B/a 3 , where N is the total 
number of links in the chain, a stands for the typical distance between the beads-monomers, and B is the second 
virial coefficient [46,47]. We consider below only the case of a polymer coil, corresponding to z > 0. In the region of 
z< 1, the perturbation theory in powers of z can be developed and for the short-range potentials one can find [46,47] 
that 

a 2 =a 2 {z) = l + k lZ + k 2 z 2 + fciwl.28, k 2 = -20.8. (67) 

One of the important problems in the physics of polymer coils, consists in the continuation of the expansion (67) 
to the region of z ~ 1 [46,47]. On the other hand, in the limit of z 3> 1, a is related to z by a simple power-law 

a^z 2 "- 1 , (68) 

where the critical index v > 1/2, can be calculated by different methods [46-48]. We propose below, using the self- 
similar renormalization, a simple way to continue (67) to the region of arbitrary z, including both known limiting 
cases and allowing to estimate v from the stability condition. The problem of this type was already mentioned above, 
in Sec. II. From the viewpoint of the applicability of the stability conditions, it is worthwhile to study a~ 2 (z) = a(z), 
re-expanding it in powers of z, so that 

a(z)w l + b lZ + b 2 z 2 + bi = -1.28, b 2 = 22.438. (69) 

The set of approximations to a(z), including the two starting terms from (69), can be written down as follows: 

a = 1, 



oi = 1 + hz, 

and the expression for the renormalized quantity al can be readily obtained: 

al = (^j-Y 1 => (^-) S1 z- Sl (z - oo), (70) 
Si — u\Z —0\ 

where the stabilizer si should be positive, if we want to reproduce in the limit of z — > oo, the correct power-low behavior 
of the a 2 (z). A different set of approximations, not including into the renormalization procedure the constant term 
from (69), has the form: 

oT = b\z, 
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a 2 = b\z + b 2 z 2 

and applying the standard procedure, we obtain 

b 2 Z l-(l+s 2 ) . 



l + Ml 



b 1 (l + s 2 ) 1 ' v l + s 2 

Demanding now, that both (70) and (71) have the same power-law behavior at z 

s 2 =si = 2(2v - 1). 



( T ^)-( 1 + s ^ b\ +s *z~ S2 (z-»oo). 

— > oo, we find that 



(71) 



Requiring now the fulfillment of the stability criteria for the two approximations (70) and (71) in the form of the 
minimal-difference condition (see Introduction), we obtain the condition on the positive stabilizer si, i.e, Si should be 
determined from the minimum of the expression A : 



U l + si 1 -&i 



(72) 



The minimum of (72) does exist and is located at the point S\ w 0.5, Correspondingly, the index v is equal to 0.625, 
in a reasonable agreement with all other theoretical and experimental estimates of this index [46-48] . 

As it was pointed out in Sec. 2, the results may also be obtained, if the self-similar renormalization is applied to the 
initial expansion (67), for the sought function a 2 (z), although it is formally divergent at z — > oo and the stabilizer 
should become negative to describe this divergence correctly. By simple substitution of the coefficients and changing 
the criteria on minimum of (72) to the maximum of the analogous expression K : 



K 



{( _^ r a +s) fc ( 2+ *)_ ( ^_ n 



(73) 



with respect to the now negative stabilizer s = 2(1 — 2v). One can see that the maximum is located at the point 
s = —0.3719, leading to the very reasonable estimate for the critical index v = 0.593. Final formulae have the following 
form: 



2/ \* 

a t {z) 



' s — k\z 



(74) 



Both formulae, (74) and (75), with s w —0.3719, can be used as an approximate "equation of state" for the polymer 
in the whole range of the parameter z, satisfying, by design, both known virial and scaling asymptotic expressions. 



VIII. CRITICAL TEMPERATURE OF THE 2D ISING MODEL FROM THE EXPANSION AROUND 

THE DIMENSION ONE 

In this Section we calculate the critical temperature T c of the two-dimensional (2D) Ising model starting from 
the approximate expression obtained within the framework of the so-called quasi-chemical approximation [42] . This 
approximation gives T c as a function of the coordination number z : 

T ^ ^ w^wr (76) 

The expression (76) correctly describes the limit of the one-dimensional (ID) Ising model, with T c = 0, and at the 
infinite-dimensionality the limit coincides with the well known Bragg- Williams result T c — z. The expansion around 
the latter limit has been widely used, although its accuracy is not too good [49]. We adopt the different approach, 
expanding (76), in powers of the parameter z — 2 = A, around its correct, D = 1 (z = 2), limit, and use the self-similar 
renormalization to continue the expansion valid at A < 1 to the region of A = 2, corresponding to the 2D Ising 
model with the quadratic lattice. The expansion of the inverse expression (76), up to the quadratic term in A, has 
the following form: 
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r-(A) ^[i + E L + Ia-^. (7-) 



We renormalize separately the logarithmic 



i(A) = [l + ^ln(l)] (78) 



and power-law contributions 

P(A) = iA-^A 2 , (79) 

separating in this way, the effects of long-range and short-range contributions to T c . 

The standard prescriptions of Sec. Ill are fully applicable to the expression (78) containing the logarithmic term. 
Two consecutive approximations to L(A) are 

MA) = 1, 



MA) = l + ^ln(i), 
with the expansion function / = A s , and the multiplier 

mi (A, S ) = l-^[ln(A) + i], 
equal, at A — 2, to zero as s — > oo. The velocity function has the following form: 

W(S ' /) = -ln2— • 
Calculating the evolution integral and taking the limit of s — > oo, we obtain 

L*(A) = A" 1 /' 112 

The expression for P(A)*, can be readily written down in the case of summation along the stable trajectory, 
corresponding to s — > oo: 

P*(A) = iAcxp(-iA). 

Recalculating T* , we obtain T*(A = 2) = 2.321. The percentage error equals 2.292% , when compared to the exact 
Onsager result T c = 2.269. It should be remembered, that the quasi-chemical approximation (76) works with the 
percentage error of 27.149%, and that the initial expansion (77) deviates from the exact result by 76.289%. Also, one 
of the best known approximate theoretical schemes, known as the Kikuchi method [42], gives the percentage error 
equal to 6.831%. 



IX. TEMPORAL ASYMPTOTES OF THE DIFFUSION EQUATION WITH RANDOM STATIONARY 

NOISE 



A. Poisson spectrum 



Consider the diffusion of particles in a medium with randomly distributed traps, whose local density a(r) is described 
by the nonnegative Poisson random field [50-52]. The local particle density n(r, t) in the presence of traps is described 
by the following equation: 

- ^n(r, t) = -V 2 n(r, t) + a(r)n(r, t), (80) 
where we put the diffusion coefficient equal to one. 
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This problem is formally equivalent to the Shroedinger equation with imaginary time and potential a(r) [50-52]. The 
eigenvalues of the Shroedinger operator corresponding to (80), Ei, are nonnegative. In addition, the density of states 
p(E), near the finite fluctuational limit of the spectrum, located at E = 0, and formed due to the rare fluctuations of 
the potential with a close to zero in large regions of space, is known explicitly: 

p(E) = f(E)cM-E- D/2 ), (81) 

where D is the dimensionality of space [53,54]. Hereafter, for the sake of simplicity, we omit the constant terms in the 
exponential for the density of states. Also, for D = 1, the pre-exponential factor f(E) ~ E~ 3 / 2 , is known [53] and, 
because of this, we consider below only the one-dimensional case. The general solution of (80) can be readily written 
down [52] in terms of the cigenfunctions and eigenvalues of the corresponding quantum-mechanical problem: 

n(r, t) = CiMr) exp(-i^), (82) 

i 

and, considering only contributions from the rare fluctuations of the potential, the mean density (n) over the entire 
volume for t — > oo can be represented in the form of the integral: 

roc 

(n{t))=n p(E)exp(-Et)dE, (83) 
Jo 

where no describes the initial homogeneous particle distribution. This integral can be evaluated by the method of 
steepest descent and the leading exponential term ~ exp(— t 1 / 3 ) can be written down [49-51]. Similar estimates were 
obtained for arbitrary D [50-52]. 

We will obtain below the higher-order contributions to the asymptotic expansion of (83) near the saddle point 
E = (^) 2 / 3 and, using the self-similar renormalization, obtain the leading corrections, as t — > oo, to the pre- 
exponential factor f(E) ~ t -1 . 

Let us represent E in the vicinity of E as E = E + e and expand the expression 

$(t,e) = ln\p(E)exp(-Et)] 
in powers of e up to the third order terms, so that 

$(t, e) w -at 1 ' 3 - A(t)e 2 + B(t)e 3 , 

a = 3 x 2- 2 / 3 , A(t) = 3 x 2- 4 / 3 t 5 / 3 , B(t) = 5 x I^IH'" 3 
and, expand the exp{B(t)e 3 } in powers of e, so that 

exp{$(t, e)} w exp{-ai 1/3 } exp{-A(t)e 2 } x [1 + B(t)e 3 + ...]. 

Now, (n(t)) can be written down as follows: 

(n(t)) - t 1/6 exp{-at^ 3 }[l + br x /% b = 0.684, (84) 

i.c, the corrections to the pre-exponential factor are obtained in the form of an expansion in inverse powers of t, valid 
as t — > oo. Our aim is to continue this expression to the region of t ~ 1. 

Apply now the self-similar renormalization to the quantity n(t) = 1 + &t -1 ^ 6 , with the two consecutive approxima- 
tions: 

n = 1, 



ni - l + 6t~ 1/6 . 

Following the standard prescriptions, the renormalized quantity n* (t) can be obtained: 



n*(t) = {i 1 /6 + _A_ } 6 S (t) ; (g5) 
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where the stabilizer 

- b 
S[t) ~ 6tV6 + 6&' 

is defined as zero of the multiplier 

mi(s,t) = 1 + % -. 

s 

For the intermediate region 1 <C t < oo, the simple expression can be written down: 

n*(t)oci"/ 6tl/6 , 

which gives the correction to the pre-exponential factor in the form of continued noninteger powers. It is worth noting, 
that already the starting terms of the asymptotic expression (85), lead to the approximation-cascade trajectory with 
zero multiplier. 



B. Gaussian spectrum. 

Consider equation of the same type that (80), with the only difference that the Poisson potential a is replaced by 
the random potential U(r) with the properties of the Gaussian "white noise", i.e. U(r) = 0, (U(r)U(r )) oc 6(r — r ) : 

- -|n(r, t) = -V 2 n(r, t) + U (r)n(r, i). (86) 

Equation of this type, but with a noise dependent both on space and time, can be easily transformed to a nonlinear 
Burgers equation, Kardar-Parisi-Zhang (KPZ) equation, and also it describes some other closely related physical 
problems [55]. Stationary random potential is not meaningless within the framework of, say, KPZ where one can 
think about the stationary, random in space, perturbations of a growing interface. It is also considered in biology as a 
model for population dynamics in the presence of a random distribution of food [55]. Equation of the type (86) with 
a random potential U (r) can be also transformed to the corresponding Shroedinger equation with imaginary time. 
Therefore, one can use the known properties of the spectrum of the corresponding quantum-mechanical problem to 
study the t — > oo behavior of the diffusion equation with stationary, randomly distributed sources and sinks. The 
fluctuational limit of the spectrum is situated now at E — > — oo, with the exponentially small density of states in its 
vicinity 

p(E)=f(E)cxp(-\E\ 2 - D ^), (87) 

being formed due to the rare fluctuations of the potential with large negative values, separated from each other by 
the distances much larger than their own sizes [53]. Such fluctuations may be again considered separately. The 
eigenvalues corresponding to the eigenfunctions localized at these fluctuations are now negative, in distinction from 
the case considered above and the mean density evolution as t — > oo can be estimated from the following integral: 

(n(t)} =n p(E) exp(| E \ t)dE. (88) 
Jo 

We again use the method of steepest descents and follow literally the same steps as in the previous Subsection. The 
saddle point does exist for 1 < D< 4 (except at D — 2, D — 4, where the situation becomes trivial) and is given by 
the expression 

|l = (JL)2/(2-D) . 

1 1 4-D 1 

The leading exponential term has the form 

<n(t)> ~ exp{a(D)t£Z}, a(D) = (|^§)(2 - f )^, 
with radically different behavior for D = 1 and D = 3 : 



19 



(n(t)) ~exp{|a(l) \t 3 }, D = 1, (89) 



(n(t)) - exp{- | o(3) | i" 1 }, £> = 3, (90) 

corresponding to an anomalously fast growth, compared to the simple exp(£), and, anomalously slow decay, compared 
to exp(— t), respectively. We think that this difference takes its origin from the principally different properties of the 
corresponding Shrocdinger operator, where it is known, that at D = 1 all states of the particle are localized, while 
at D = 3, generally speaking, both localized and delocalized states are present [53]. These basic theorems, when 
applied to the case of diffusion, explain why at ID , the random distribution of sources and sinks, causes an explosive 
instability of the density fluctuations, while in 3D, the disorder can cause only longer decay times for the density 
fluctuations. Of course, the instability can be cured by nonlinear interactions, that should be now taken into account. 

At D = 1, where f{E) ~ E [54], applying the procedure already discussed above, we obtain the expansion in the 
vicinity of the saddle point: 



(n(t)) oc cxp{| o(l) | t 3 } x i 5/2 [l + bt- 3/2 + ...], 




and the renormalized expression can be readily written down, using the same definition for n{t) as above: 

At D = 3, a little bit different situation occurs, since 

$(t, e) w - | o(3) | t- 1 + A(t)e 2 - B(t)e 3 , 

A(t)=t 3 , B(t)=2t 5 , 

and, in order to guarantee the convergence of the integrals, it is reasonable to expand exp($(t, e)) as follows: 

cxp{$(t,e)} w cxp{- | o(3) | t _1 }exp{-B(t)e 3 } x [1 + A(t)e 2 + ...], 
and for the (n(t)) we obtain 

(n(t)) - cxp{- | a(3) | i" 1 }*- 5 / 3 ^ + rf- 4 / 3 + ...], c w 0.296. 
The renormalized expression for the pre-exponential factor has the following form: 

We have demonstrated in this Section, that the self-similar renormalization can be applied to the dynamical problems 
as well, generating the expressions for the pre-exponential factors in the form of continued noninteger powers. 

X. CONCLUSION 

We suggested here a new variant of the self-similar approximation theory permitting an easy and accurate summation 
of divergent series. The method is based on a power-law algebraic transformation leading to an effective increase of 
the order of perturbative terms. The powers of this transformation play the role of control functions governing the 
convergence of renormalized series. These control functions arc defined by the principle of maximal stability, i.e., from 
the minimization of mapping multipliers. Such stabilizing control functions may be called stabilizers. 

Another important novelty of the method is the multiple self-similar renormalization converting all series into 
closed self-similar expressions. This multiple and complete renormalization is called self-similar bootstrap. The 
resulting effective sum of a divergent series can be presented through analytical expressions containing exponentials 
and fractions, rational or irrational. In particular cases, this can be only exponentials or only fractions, depending on 
the behavior of control functions which dictate the resulting form. Because of much larger variety of such resulting 
forms, the method allows to present the answers in relatively simple analytical expressions having at the same time 



20 



quite high accuracy. The use of several types of functions, such as exponentials and various fractions, distinguishes 
this method from, e.g., the Pade approximants involving solely rational functions. 

In order to prove that the suggested method really gives quite simple and accurate expressions for the effective 
sums of divergent series, we , first of all, considered several toy models cartooning the generating functionals in field 
theory or partition functions in statistical physics. We illustrated by these examples that the method works well in 
different situations, for single-well and for double-well models, for weak and strong coupling. 

To stress the generality of the method, we applied it to several problems of statistical physics of quite different nature: 
to constructing the equation of state, to calculating the critical temperature, and to finding the time asymptotics 
for stochastic dynamical processes. We hope that these various and very different applications demonstrate well the 
validity of the method. 
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